Vortex phase diagram for mesoscopic superconducting disks 
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Solving numerically the 3D non linear Ginzburg-Landau 
(GL) equations, we study equilibrium and nonequilibrium 
phase transitions between different superconducting states of 
mesoscopic disks which are thinner than the coherence length 
and the penetration depth. We have found a smooth tran- 
sition from a multi-vortex superconducting state to a giant 
vortex state with increasing both the disk thickness and the 
magnetic field. A vortex phase diagram is obtained which 
shows, as function of the magnetic field, a re-entrant behav- 
ior between the multi-vortex and the giant vortex state. 

PACS number(s): 74.24.Ha, 74.60.Ec, 73.20.Dx 



Recently, mesoscopic superconductivity has attracted 
much attention in view of phase transitions in confined 
systems with sizes comparable to the coherence (£) and 
penetration (A) lengths. While the type of bulk supercon- 
ductors is only determined by the value of the Ginzburg- 
Landau parameter k = A/£, the experimental observa- 
tions and the numerical simulations of magne- 
tization of mesoscopic thin disks have shown that the 
type and the order of those transitions between differ- 
ent superconducting states and between the supercon- 
ducting and the normal state depends crucially on the 
disk radius R and the thickness d. With increasing the 
disk radius the second-order reversible phase transition 
observed for small disk radii are replaced by first-order 
transitions with a jump in the magnetization. In pre- 
vious theoretical investigations of superconductivity in 
such mesoscopic disks only the giant vortex states 
with fixed total angular momentum L were considered 
with an axially symmetric order parameter. It is well 
known [jjj that for type-II superconductors (k > l/\/2), 
the triangular Abrikosov vortex lattice is energetically 
favorable in the range H c i < H < _ff c2 . Since the ef- 
fective London penetration depth A = X 2 /d increases 
considerably in thin disks with d -C A one would ex- 
pect the appearance of the Abrikosov multi-vortex state 
even in disks made from a material with k < 1/ \/2, like 
e.g. the Al disks studied in Refs. |^,|3|]. By analogy with 
classical particles confined by an external potential 
the structure of a finite number of vortices should dif- 
fer from a simple triangular arrangement and allow for 
different metastable states. Using the London approxi- 
mation Fetter || calculated the critical field H c \ for flux 
penetration into a disk. For a superconducting cylinder 
the multi-vortex clusters, containing up to four vortices, 
were simulated by Bobel 0. Using the method of im- 
ages and the London approximation, Buzdin and Brison 



H have considered vortex structures in small R <C A 
disks and found a classical particle ringlike arrangement 
H of vortices. In the present Letter we study the tran- 
sition from the giant vortex state to this multi-vortex 
configuration for thin superconducting disks within the 
nonlinear Ginzburg-Landau (GL) theory. 

We consider a superconducting disk immersed in an 
insulator media with a perpendicular uniform magnetic 
field H . For thin disks (d <C £, A) we found (||| that 
it is allowed to average the GL equations over the disk 
thickness. Using dimensionless variables and the London 
gauge divA = for the vector potential A, we write the 
system of GL equations in the following form 
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with the boundary condition (— iV 2 £> — j4) v E , | r =7? = 0. 
Here the distance is measured in units of the coherence 
length £, the vector potential in c?i/2e£, and the mag- 
netic field in H C 2 = ch/2e^ 2 = k\/2H c . The disk is 
placed in the plane (x, y), the external magnetic field is 
directed along the z-axis, the indices 2D, 3D refer to two- 
dimensional and three-dimensional operators, and jiD is 
the density of superconducting current. To solve the sys- 
tem of Eqs.(^jj^) we apply a finite-difference representa- 
tion of the order parameter and the vector potential on a 
uniform cartesian space grid (x. y) and use the link vari- 
able approach || . In our simulations we took typically a 
grid spacing of 0.15£. To find the steady-state solution 
of the GL equations we add to the LHS of Eqs. (Q) and 
(||) the time derivatives of the order parameter and the 
vector potential, respectively, and use an iteration proce- 
dure based on the Gauss-Seidel technique to find *. The 
vector potential is obtained with the fast Fourier trans- 
form technique. For this purpose we set the condition for 
the vector potential A\ x \—n s i„i_^ s = Hq(x, —y)/2 at the 
boundary of a larger space grid (R s = 4i?), where the 
vector potential created by the superconducting currents 
is much less than the external vector potential. 

The giant vortex state is characterized by the total an- 
gular momentum L through ^ = ?p(p)exp(iL(j)), where p, 
tf> are the cylindrical coordinates. An arbitrary supercon- 
ducting state is generally a mixture of different angular 
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harmonics. Nevertheless, we can introduce an analog to 
the total angular momentum which is still a good quan- 
tum number. Choosing circular loops at the periphery 
of the disk we find that the effective angular momentum 
L — A0/2-7T does not depend on the loop radii pi when 
it is in some range pi = (0.8 -j- 1)R. This allows us to 
characterize unambiguously the different superconduct- 
ing states. Note, that the effective angular momentum 
is in fact nothing else then the number of vortices in the 
disk. 

To find the different vortex configurations, which in- 
clude the metastable states, we search for the steady- 
state solutions of Eqs.(Q-||) starting from different ini- 
tial conditions which were generated randomly. Once 
we obtained these configurations for some value of the 
magnetic field we increase/decrease slowly the magnetic 
field as long as the number of vortices remains un- 
changed. Comparing the dimensionless Gibbs free en- 
ergies F = V 1 J(2(A - A )j2d - |*| 4 )dr, where inte- 
gration is performed over the disk volume V, and Aq 
is the vector potential of the external uniform magnetic 
field, we found the ground state. If the system has suffi- 
cient time to equilibrate, the system will be in the ground 
state and we obtain phase transitions between different 
ground states. If the latter condition is not satisfied the 
system can remain in a metastable state until this state 
disappears or becomes unstable with respect to small per- 
turbations in the order parameter and/or the magnetic 
field. Such nonequilibrium phase transitions lead to hys- 
teresis in the magnetization. The barriers separating the 
metastable and the ground state are discussed in Refs. 
0,0. The system remains in the metastable state if 
the barrier height is much larger than the temparature. 
Instead of finding the barrier heights, which is a more 
cumbersome problem, we search for the critical magnetic 
fields corresponding to the disappearance of such barri- 
ers. Here, we started from the Meissner (normal) state 
and increased (decreased) slowly the magnetic field (typ- 
ically with steps of O.Oli/^) from the zero (nucleation) 
magnetic field. 

The free energies of the different vortex configurations 
are shown in Fig. [j] for zero disk thickness and for two 
disk radii: (a) R = 4£, and (b) R = 4.8£. For R = 4£, the 
vortex configuration can consist up to six vortices which 
are arranged on the edge of an ideal polygon. The pen- 
tagon and hexagon vortex clusters are always metastable 
states for R = 4£. With increasing disk radius the al- 
lowed number of vortices increases leading to the appear- 
ance of polygons with a vortex inside it when L > 7. 
We found that these structures have a larger energy as 
compared to the ideal polygons. This result differs from 
the London approximation of Ref. Q where the vortices 
were treated as an ensemble of interacting classical point 
pseudoparticles in which the closed packed structures are 
more preferable even for L > 6. The present work clearly 
shows the limited validity of the approach of Ref. || for 
the case of superconducting disks. Another unexpected 



feature is observed when the magnetic field is further in- 
creased. While the model of classical particles predicts a 
decrease in the intervortex distance followed by the ap- 
pearance of a new vortex, our simulations show, as a rule, 
a gradual transition from a multi-vortex state to a giant 
vortex state (Fig. ||). In principle, the latter may be a 
metastable state. The positions of these transitions are 
indicated by circles in Fig. |l|. For a larger disk radius 
R = 4.8£, we also observe transitions between different 
multi-vortex states. Performing a Fourier analysis of the 
order parameter we find that a multi-vortex state corre- 
sponding to an ideal polygon presents a mixture of har- 
monics ^ pa ipo(p) + Tlik' l ^k{p)e.xp{ikL(j)) with a rather 
small contribution of higher (k > 2) harmonics. This al- 
lows us to find approximately the free energy of different 
vortex configurations in thin disks i? < A in which we 
can neglect the distortion of the magnetic field. For this 
purpose we take the order parameter as a superposition 
of only two states ~ C^ 2 (o(p) + C\l 2 C,L{p)exp{iL(j)), 
where (o , Cl arc the cigenfunction of the linearised first 
GL equation (|l|) for different angular momenta P,p^[. 
Substituting this expansion into Eq. (|I|) we obtain the fol- 
lowing set of non-linear equations for the coefficients Co 
and Cl 

XqCq = anC 2 + a 12 C C L , (4) 
AlCl = ouCqCl + C122CI, 

with an =< Co ICo >> a i2 = 2 < Cold >, «22 =< 
ClICl >j where Ao, Xl are the eigenvalues of the lin- 
earised GL equation ||, and < /o|/i > refers to the 
matrix element V^ 1 Jfofidr. Besides the two triv- 
ial solutions Cl = 0, Cq = Ao/aii and Co = 0, 
Cl = Xl/o22 which correspond to the Meissner state 
and the giant vortex state, respectively, it is possible 
to have another solution with Co = (Ao£t22 — Xlo^/D, 
Cl = (Al<Zii— Aoai2)/A D = 0,110,22— a\ 2 corresponding 
to a multi-vortex state. The free energy of the ground 
state obtained using this approach, is shown in Fig. [I] by 
the dashed curve, which is in good agreement with the re- 
sults of our simulations (full curves). In order to discuss 
the multi-vortex <-> giant vortex transition, we anylizcd 
the stability of the obtained solutions with respect to 
small perturbations of the coefficients Co, Cl (see Ref. 
||). The stability conditions for the multi-vortex state 
and the giant vortex state are Xl < A* = AoG,22/ai2 and 
Xl > A*, respectively. Consequently for fixed L, there is 
an unique solution with a reversible transition from the 
multi-vortex state to the giant vortex, which occurs with 
increasing magnetic field when = A*. At this critical 
point the multi-vortex state coincides with the giant vor- 
tex state and consequently there is no jump in the mag- 
netization. But the derivatives of the coefficients Co, Cl 
are discontinuous and correspondly we expect a discon- 
tinuty in the first derivative of the magnetization. This 
is also confirmed by our numerical simulations where in 
Fig. |3| the region around the L = 3 multi-vortex to giant 
vortex transition is shown. 
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The magnetization of the disk M — J (H — 
H )df/4:TTVH C 2 with R = 4.8£ is shown in Fig. |J for 
increasing (a) and decreasing (b) magnetic field. Jumps 
in the magnetization correspond to transitions between 
states with different number of vortices at the magnetic 
fields where the states cease to exist (see Fig. 1). With 
increasiing disk thickness the demagnetization effect in- 
creases and as a result the transition points shift to larger 
magnetic fields, but the number of jumps remains the 
same. The average magnetic field in the disk can be 
estimated as (H) w H — AttM and the relative magneti- 
zations M/d as a function of (H) is 'almost' an universal 
curve for all disk thicknesses (Fig. ||(c)). For increasing 
magnetic field and disk thicknesses d — 0.2, 0.4£ all tran- 
sitions occur between giant vortex states. For the thinner 
disk d — 0.1£ we observe the following sequence of transi- 
tions -> f -> 2g — > 3 m — > 3g — > 4 m -> 4g -> 5g ~> 6g . . . , 

where the lower index (g,m) corresponds to giant and 
multi-vortex states, respectively. For decreasing mag- 
netic field multi-vortex states appear with a larger num- 
ber of vortices. The sequence of transitions between 
multi-vortex states start from L = 7 and L = 8 for 
d = 0.2£ and d = 0.1£, respectively. For the thicker disk 
(d = 0.4£), a multi-vortex state appears only for L = 3, 2 
just before the expulsion of a vortex from the disk. As is 
evident from Fig. ^(b), the appearance of positive mag- 
netization observed experimentally can also be explained 
within the 3D GL approach without the consideration 
of pinning effects. But the lowest energy state, i.e. the 
equilibrium state, has always a negative magnetization. 

To distinguish quantitatively the giant vortex state 
with the multi-vortex states with the same number of 
vortices we consider the value of the order parameter 
l^j 2 in the center of the disk which is nothing else than 
the density of Cooper pairs. We find that this param- 
eter, which is zero for a giant vortex state, goes almost 
linearly to zero when the magnetic field approaches some 
critical value (see Fig. || (a), the thick curve). There- 
fore, the magnetic field obtained by linearly interpolat- 
ing 1^(0, 0)| 2 to zero defines the transition from a multi- 
vortex state to a giant vortex state. As mentioned above, 
our simulations show that this transition is reversible and 
exhibits a discontinuous derivative in the magnetization. 
Having the free energies of different vortex configurations 
we are able to construct an equilibrium vortex phase di- 
agram which is shown in Fig. |^ for two disk radii R = 4£ 
and R = 4.8£, respectively. The solid curves separate 
the regions with different number of vortices and the 
dashed curves show the boundaries between the multi- 
vortex and the giant vortex states. For L = 1 the sin- 
gle vortex state and the giant vortex state are identi- 
cal. The shaded regions correspond to the multi-vortex 
states. The superconducting to normal transition occurs 
for H/H C 2 w f.9 which is outside the plotted region. 
Notice that the multi-vortex area in the phase diagram 
reduces in size with incresing disk thickness and it dis- 
appears in the limit of thick disks where only the giant 
vortex state survives. Thus for type 1 superconductors 



the multi-vortex state is favored with decreasing disk 
thickness and increasing disk radius. This behaviour can 
be understood as follows: with increasing radius the en- 
ergy difference between different L-states decreases and 
consequently it becomes possible to build a lower energy 
multi-vortex state out of a linear combination of giant 
vortex states. For decreasing radius this is more difficult 
to do and there exists a critical radius below which no 
multi-vortex states have the lowest energy. To observe 
the multi-vortex <-> giant vortex transition one should 
investigate the derivative of the magnetization or the vor- 
tex configuration itself which can be done by using e.g. 
a magnetic force microscope. 
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FIG. 1. The free energy of configurations with differ- 
ent number of vortices L for two disk radii R = 4£ (a) and 
R = 4.8£ (b) in the case of zero disk thickness and k = 0.28. 
The open circles indicate the transition from a multi-vortex to 
a giant vortex state (right side of circle). The dotted curves 
correspond to the free energy of metastable configurations 
with a vortex inside a poligon for L = 7, 8 (b). The dashed 
curves are the results of our approximate analytical calcula- 
tions. The insets show the possible vortex configurations. 

FIG. 2. Contourplot of the magnetic field distribution in 
the disk plane (z = 0) for the case of a three vortex state 
and for different applied magnetic fields Ho = 0.525_ff C 2 (a), 
H = 0.65H c2 (b), H = 0.75H c2 (c), and H = 0.8H c2 (d). 
We took R = 4£ and d — 0.5£k 2 . For the lowest magnetic 
field (a) the three-vortex configuration is metastable and a 
small decrease in the magnetic field leads to a transition to 
the two vortex state due to the expulsion of a vortex through 
the disk boundary indicated by the arrow in (a). 

FIG. 3. The free energy (I), the square of the order pa- 
rameter in the disk center (2), the magnetization (3), and its 
first derivative (4) for a disk with R = 4£, d — 0.5£k 2 which is 
in the vortex state with L — 3. The dashed vertical line shows 
the transition from the multi-vortex state to the giant vortex 
state. In (a) the solid circles indicate the points at which the 
equilibrium phase transitions 2 — > 3 and 3 — * 4 occur. 

FIG. 4. Magnetization of a disk in increasing (a,c) and 
decreasing (b) magnetic fields for R — 4.8£, n = 0.28 and 
different disk radii. 



FIG. 5. The vortex phase diagram for two different disk 
radii R = 4£ (a) and R = 4.8£ (b). The shaded area corre- 
sponds to the multi-vortex state 



4 




I I I I I I I I I I I I I I I I I I I I l_ 

-4 -2 2 4 -4 -2 2 4 

x/5 x/5 



Fig.2 




0.0 0.5 1.0 1.5 2.0 

H/H c2 

Fig.4 




Fig. 5 



